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ABSTRACT 

We use the probability distribution function (PDF) of the Lyman-a forest flux at 
z = 2 — 3, measured from high-resolution UVES/VLT data, and hydrodynamical 
simulations to obtain constraints on cosmological parameters and the thermal state 
of the intergalactic medium (IGM) at z ~ 2 — 3. The observed flux PDF at z = 3 
alone results in constraints on cosmological parameters in good agreement with those 
obtained from the WMAP data, albeit with about a factor two larger errors. The 
observed flux PDF is best fit with simulations with a matter fluctuation amplitude 
of as — 0.8 — 0.85 ± 0.07 and an inverted IGM temperature-density relation (7 ~ 
0.5 — 0.75), consistent with our previous results obtained using a simpler analysis. 
These results appear to be robust to uncertainties in the quasar (QSO) continuum 
placement. We further discuss constraints obtained by a combined analysis of the 
high-resolution flux PDF and the power spectrum measured from the Sloan Digital 
Sky Survey (SDSS) Lyman-a forest data. The joint analysis confirms the suggestion 
of an inverted temperature-density relation, but prefers somewhat higher values (ag ~ 
0.9) of the matter fluctuation amplitude than the WMAP data and the best fit to 
the flux PDF alone. The joint analysis of the flux PDF and power spectrum (as 
well as an analysis of the power spectrum data alone) prefers rather large values for 
the temperature of the IGM, perhaps suggesting that we have identified a not yet 
accounted for systematic error in the SDSS flux power spectrum data or that the 
standard model describing the thermal state of the IGM at z ~ 2 — 3 is incomplete. 
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1 INTRODUCTION 

The Lyman-a forest is an important cosmological observable 
that probes matter density fluctuations in the IGM over a 
unique range of redshifts, scales and environments. Many 
attempts have been made to measure physical properties of 
the IGM using Lyman-a forest data. The two most common 
approaches are either based on decomposing the information 
encoded in the transmitted flux via Voigt profile fitting or 
treating the flux as a continuous field with directly mea- 



surable statistical properties (e.g. Rauch et alj 1997: 



1 19981 : iTheuns et al] Il998l : I Croft et alj |2002t lMeiksinll2007l ) 



Ranch 



In the second approach, measurement of the zero, one, two- 
point or three-point probability distribution functions (i.e. 
the mean flux level, the flux PDF, the flux power and bispec- 
trum) enable a variety of physical properties to be explored. 
The mean flux level for example, is sensitive to the am- 



plitu d e of the meta-gala ctic UV background (|Tvtler et alj 
12004 iBolton et all 120051 ) while the flux PDF is sensitive 
to the thermal evol ution of the IGM l|McQuinn et alj|2009l : 
IBolton et al . 2009b). The flux power spectrum has been used 
to constrain cosmological par ameters and the behaviour of 
dark matter at small scales (jViel et alj 120041 : ISeliak et al.l 
20061: IViel et alJfeOQST ) and the flux bispectrum can be used 
to search for signatures of n on-gaussianities in the matter 
distribution l|Viel et al.ll2009f ). Ideally, a given IGM model 
described by a set of cosmological and astrophysical param- 
eters should agree with all these statistics including the re- 
sults from Voigt profile decomposition at the same time. In 
practice, the interpretation of the data is complex and is 
heavily dependent on numerical simulations that incorpo- 
rate the relevant physical ingredients, but have a limited 
dynamic range. 
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The data used for these investigations consist mainly 
of two kinds of sets of QSO spectra: the SDSS low- 
resolution, low signal-to-noise sample and UVES/VLT or 
HIRES/KECK samples of high-resolution spectra. The 
characteristics of the low and high-resolution data sets are 
very different (the number of SDSS spectra is about a fac- 
tor ~ 200 larger than that of high-resolution samples, but 
the latter probes smaller scales due to the higher spectral 
resolution). Measurements based on Lyman-c* forest data 
have reached a level of accuracy where an understanding of 
systematic uncertainties at the percent level or below (the 
magnitude of statistical errors associated with the SDSS 
sample) has become important. In this Letter, we will re- 
visit the flux PDF, which has been investi gated previously 
by several authors either on its own (e.g . 
200d; I Jena et al.|l2005l; IBecker et al.ll2007l ; 
Kim et al.ll2007l ; iBolton et alJl2o6i T or jointly with the flux 



McDonald et all 
Lidz et al.ll2006l ; 



power spectrum (iMeiksin et al.l l200ll : IZaroubi et al.l 120061 : 
iDesiacaues et al.ll2007l '). We shall focus on the fl ux PDF of 
the U VES/VLT data as recently measured by iKim et all 
1)2007! ) (hereafter K07); the systematic and statistical er- 
rors for this sampl e have been addressed in detail. Recently, 
bot h IBolton et all d200Sh (from the K07 UVES/VLT data) 
and IBecker et alJ (|2007f ) (from independent HIRES/KECK 
spectra) have found evidence for a density-temperature re- 
lation of the IGM which appears inverted if approximated 
as a power-law (7 < 1 for T = To(l + ( S) 7 " 1 ). Here, we wil l 
improve on the analysis performed in IBolton et al.l (I2008T ) 
(hereafter B08) and check its robustness by fully exploring 
the cosmological and astrophysical parameter space. In ad- 
dition, we briefly discuss a joint analysis of the flux PDF 
and SDSS flux power spectrum, and the possible implica- 
tions for constraints on cosmological parameters describing 
the linear matter power spectrum and the thermal history 
of the IGM. 



2 METHOD 

We use simulations performed with the parallel h ydrody- 
namical (TreeSPH) code GADGET-2 (|Springelll2005h to cal- 
culate the flux statistics for models with a wide range of 
cosmological and astrophysical parameters by expanding 
around a reference model. For the reference model we choose 
here the 20-256 simulation of B08. We refer the reader to 
this paper for further detai ls, including resolution and box 
size convergence tests (see IBolton fc Beckerl (|2009h for re- 
cent convergence tests on SPH simulations). We will com- 
pare these simulations to improved measurements of the 
PDF made by K07 in three redshift bins at (z) = 2.07, 
(z) — 2.52 and (2) = 2.94 based on a set of 18 high reso- 
lution (R ~ 45 000), high signal-to-noise (S/N > 30 - 50) 
VLT/UVES spectra. Further details regarding the observa- 
tional data and its reduction, with particular emphasis on 
metal removal and continuum fitting errors, may be found 
in K07. In all instances the mock QSO spectra have been 
processed to have the same instrumental properties as the 
observed data: i.e. the same signal-to-noise, resolution and 
pixel size. 

We explore the following cosmological and astrophys- 
ical parameters: Q m , n s , Ho, as for the cosmological part 
3) and >y A ' s (z = 3) for the IGM thermal 



Par am. 


(20,256)°' 1 °' ^ 


(20,256)° °'^ 


(20,256)° 1 




0.81 ±0.07 (0.80) 


0.85 ± 0.07 


0.86 ±0.06 


n s 


0.96 ± 0.03 (0.95) 


0.96 ± 0.03 


0.96 ±0.03 




0.23 ±0.07 (0.19) 


0.22 ± 0.06 


0.22 ±0.06 


H 


82 ± 7 (80) 


84 ±6 


85 ±8 


T 


19 ±6 (15) 


24 ±8 


26 ±7 


qnsl 


n fi + 1 4 (— H) 


1.3 ± 1.2 


1.5 ±1.0 


A 


0.75 ±0.21 (0.72) 


0.51 ± 0.13 


0.51 ±0.13 




-1.0 ± 1.0 (-1.7) 


-1.6 ±0.9 


-1.3 ± 1.0 


T eff 


0.312 ±0.012 (0.312) 


0.321 ±0.010 


0.324 ±0.010 


T eff 


3.17 ±0.18 (3.20) 


3.16 ± 0.13 


3.25 ±0.14 


fc X 100 


± 1 (0) 


-0.8 ±0.4 


-1 ± 0.4 


X 2 /d.o.f. 


35.2/36 


45.6/48 


64/54 



Table 1. Marginalised cosmological and astrophysical parameters 
derived from fitting the flux PDF at z = 2.07, 2.52, 2.94 in the 
flux ranges F=[0.1-0.8], F=[0.1-0.9] and F=[0-1] (left, middle and 
right columns respectively): To is measured in units of 10 3 K, Ho 
in km/s/Mpc. The probabilities of having a \ 2 larger than the 
obtained values are 50, 57 and 17%, respectively. The values in 
parentheses are best-fit values for (20,256)l 01 -°' 8 ' . 



and T A,S (« 



history, where A and S indicate the amplitude and slope 
for the temperature and 7 relations normalised at z = 3 
(y = A[(l + z)/4] s ). The amplitude and slope of the effective 
optical depth evolution, t c s = — (F), are varied assuming a 
power-law evolution with redshift in order to conservatively 
span the observed range suggested by high-resolution and 
low-resolution data sets. We furthermore varied the reion- 
ization redshift (z IC — 9 in our reference model) but found 
this had no impact on the flux PDF at z < 3 (although 
the differences in "Jeans smoothing" will be important at 
redshifts close to z re , e.g. IPawlik et al. (12009 ] ')'). For the ef - 
fect on the flux power we refer to iMcDonald et all (I2005T ) 
and IViel fc Haehneltl l|2006Tl . We also consider the effect of 
a misplaced continuum level by adding an extra parameter 
f c (the flux following a continuum correction is assumed to 
be Fx (1 + /„)). 

We compute derivatives of the flux statistics from the 
20-256 model at second order using between two and four 
simulations for each cosmological and astrophysical param- 
eter. For the thermal history we explore a wide range of 
possible To and 7 values by extending the original grid of 
simulations presented in B08. The 20-256 model has 7 ~ 1.3 
below 2 = 3 and the temperature at mean density in the 
three PDF redshift bins are T = 14.8, 17.6, 20.8 x 10 3 K. 
This model was shown to be a poor fit to the K07 data in the 
simple analysis performed by B08. Here we will calculate the 
X 2 of our models varying all the parameters that affect the 
flux PDF and not just the effective optical depth, enabling us 
to expand around this mo del. This simple Taylor-e xpansion 
method was introduced in lViel fc Haehneltl l|2006l ) in order 
to explore constraints for the SDSS flux power spectrum. It 
has the advantage of enabling the exploration of the param- 
eter space close to the reference model with an accurate set 
of hydrodynamical simulations. However, the full parameter 
space cann ot be probed in this wa y with the same high ac- 
curacy (see IMcDonald et al.ll2005l for a different approach). 
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Figure 1. Best fit to the flux PDF for model (20-256) 1 " ' 8 
(continuous blue) and (20-256) 0- 1 (dashed green) in the three 
redshift bins at (z) = (2.07,2.52,2.94). 

3 RESULTS FOR THE FLUX PDF 

We obtain the best fit to the observed flux PDF for three flux 
intervals F = [0.1-0.8], F = [0.1-0.9], F = [0-1]. Different 
flux levels are subject to different systematic effects, such 
as the presence of noise and strong absorption systems at 
F ~ and the effect of continuum fitting errors at F ~ 1 
(see K07 for details). Since the PDF error bars are correlated 
we expect these systematic errors to nevertheless impact on 
the PDF over the full flux range. The level of consistency 
between the fits to these three flux intervals should indicate 
to what extent these systematic errors may or may not affect 
the results. 

For the flux range F = [0.1 — 0.8] we have a total of 
45 data points to fit and a set of 9 free parameters that 
will be varied in the Markov Chain Monte Carlo routines. 
We use the following priors on the effective optical depth, 
rjj = 0.36 ± 0.11 and r c s tf = 3.65 ± 0.21, based on the ob- 
servational results obtained by K07. Note, however, that the 
final results are affected very little by the choice of these pri- 
ors; the constraints on the effective optical depth amplitude 
at z = 3 are in fact much tighter than these priors assume. 
The results of this analysis are summarised in Table [T] and 
Figs. [T] and [2] We obtain a very good fit to the flux PDF 
for the flux range F = [0.1-0.8] (reduced % 2 = 0.98, a 
X 2 larger than this has 50 per cent probability). With the 
data points at F < 0.1 we obtain slightly larger values for 
as and Tb, but the results are in agreement at the la level. 
Adding the flux range at F > 0.9 results in a poor fit unless 
the error bars on the last two data points are increased by a 
factor of four. Note that the covariance properties of these 
data points are strongly influenced by the choice of the con- 
tinuum level. Increasing the error bars by this factor would 
account for a misplacement of the continuum level by a few 
percent. We also find evidence (2<r) that the data prefer a 
non-zero continuum offset, / c , when we add the regions at 




Figure 2. 2D marginalized and mean likelihood contours for 
cosmological and astrophysical parameters for flux PDF model 
(20,256)°' 1— 0,8 (continuous curves and filled contours). Marginal- 
ized likelihoods for the model (20,256) <)_1 are also shown in green. 



low and high transmissivity, but note that realistic errors in 
the assumed continuum level should depend on the flux and 
noise level and are expected to vary along the spectrum. 

In order to further explore the sensitivity to the last two 
data points (F = 0.95, 1) in each redshift bin, we performed 
an additional analysis by combining the two data points for 
the highest flux levels into one point at F = 0.975. We then 
recomputed the data covariance matrix and the PDF deriva- 
tives (without multiplying the covariance values of this data 
point by four). In this instance the results are consistent 
with those for the F = [0 - 0.9] and F = [0 - 1] flux ranges 
range to within la. In this case we obtain \ 2 = 60.6 for 
51 d.o.f., which indicates a reasonable fit (the probability 
for a value larger than this is 17% ). We therefore conclude 
from the results in Table [1] that our findings are robust to 
continuum fitting uncertainties and that the impact of con- 
tinuum uncertainties is mainly restricted to the flux range 
F = 0.975 - 1.025. 

The effective optical depth is constrained very well by 
the data with a best-fit value of T e «{z = 3) = 0.31 ± 
0.01, consistent with observational measurements from high- 
resolution spectra. The constraints on the temperature den- 
sity relation at z = 3 are (T A , 7 A ) = (19 ± 6 , 0.75 ± 0.21), 
in agreement with the findings of B08, and no significant 
evolution in the equation of state be l ow z < 3 is inf e rred ( in 
agreement with ISchave et all (|2000h ; iRicotti et ail (j2000h V 
The PDF alone provides interesting constraints on cosmo- 
logical parameters describing the evolution of the linear 
power spectrum; as and n s are constrained to be in the 
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range erg = 0.8 - 0.85 ± 0.07 and n s = 0.96 ± 0.03. The de- 
rived cosmological parameters are in good agreement with 
the results of other large-scale structure probes such as 
WMAP and weak lensing data, albeit with about a factor 
two larger errors t han t hose from the WMAP data (e.g. 
iLeseoureues et all (|2007n ). 

Our results corroborate the suggestion of an inverted 
temperature-density relation at z = 3. As for the analysis 
of the full flux range the statistical significance of the data 
favouring an inverted temperature-density relation 7 < 1 is 
about 3<r at 2 ~ 3. At 2 < 3 the data is consistent with 
an isothermal (7 ~ 1) temperature-density relation The re- 
sults regarding the thermal state of the IGM do not change 
significantly if we omit the flux range F > 0.9. If we dis- 
card both the flux ranges at low and high emissivity and 
consider only the flux range F — [0.1 — 0.8], there is still ev- 
idence for an inverted T — p relation, but at a reduced level 
of significance (1 — 1.5a C.L.). The likelihood contours in 
Fig. [2] indicate that a value of 7 ~ 1. 3 suggested rece n tly by 
the Hen reionization simulations of iMcQuinn et all (|2009l ) 
is between 2 and 3 — 3.5<r discrepant with the marginalized 
value we obtain when fitting the flux range F — [0.1 — 0.8] 
and the full flux range, respectively. 



Paxgrn . 


1 


B ,O»*0 
1 


R n(^ olt,2 J 1 ion 9^R~i0.1 — 0.8 


CT8 


0.85 ±0.05 


0.86 ±0.04 


0.90 ±0.02 


n B 


0.93 ±0.03 


0.96 ±0.02 


0.95 ± 0.02 




0.25 ±0.04 


0.26 ±0.03 


0.25 ± 0.03 


Ho 


78 ±7 


77 ±7 


80 ±5 


To 


38 ±7 


42 ± 6 


26 ±4 


ri-S 1 

°„ 


-0.6 ± 1.3 


-0.3 ± 1.1 


1.4 ±0.5 


rps2 


-2.3 ± 1.3 


-3.9 ± 1.3 


—3.1 ± 1.5 


°A 


0.63 ±0.50 


0.79 ±0.51 


0.50 ±0.20 


^1 


-0.7 ±2.6 


0.4 ±2 


-2.1 ± 1.6 




-1.2 ±2.1 


-1.4 ±1.5 


-1.4 ± 1.6 


T A 
'elf 


0.326 ± 0.028 


0.322 ± 0.028 


0.320 ± 0.007 


T eff 


3.19 ±0.25 


3.25 ±0.23 


3.12 ± 0.10 


2re 


11.9 ± 3.8 


9.1 ± 2.7 


14.3 ±3.6 


X 2 /d.o.f. 


78.9/85 


139.7/121 


136/130 


Table 2. 


Cosmological and 


astrophysical 


parameters derived 



from the B2i model (7 ~ 1) for the flux power spectrum: (lowz) 
flux PS fitted in the range z = [2.2 - 3.6]; (allz) flux PS fit- 
ted in the range 2 = [2.2 — 4.2]. The probability of having a \ 2 
value larger than this for model B2^ ol " z ' is 64 %, while for model 

B2j a!iz ' it is 12 %. The constraints for a joint analysis of flux PS 
and PDF (F = 0.1 — 0.8) are shown in the 3 rd column (prob. is 
35%). 



4 ADDING THE FLUX POWER SPECTRUM 

In this section we first revisit constraints from the SDSS 
flux power spectrum (PS) alone before proceeding to com- 
bine this data set with the flux PDF of the UVES/VLT data 
for a joint analysis. The SDSS flux power spectrum is based 
on 3035 QSO spectra with low resolution and low signal- 
to-noise, spanning the redshift range z = 2.2 — 4.2 (mea- 
surements are made at 11 wavenumbers in 12 redshift bins). 
Dealing with the systematic uncertainties of low resolution 
and low signal-to-noise QSO spectr a and extracting the flux 
power is a difficult task. We refer to lMcDonald et all (|2005l ) 
for a comprehensive study of the removal of continuum fluc- 
tuations, metal line contamination, damped Lyman-a sys- 
tems and dealing with the resolution of the spectrograph 
and noise level in each of the redshift bins. All these effects 
need to be properly taken into account, as a poor treatment 
would impact the obtained flux power in a non-trivial way. 
In the following, we will use the flux power provided by the 
SDSS collaboration, introducing "nuisance parameters" for 
the resolution and nois e in each redshift bin as suggested by 
iMcDonald et al l l|2005l ). and implicitly assuming that all the 
contaminants above have been either removed or properly 
modelled. 

We compute the constrai nts from the SDS S flux power 
spectrum in a similar way to IViel fc Haehn elt (2006) with 
the notable difference that we calculate the predicted flux 
statistics by expanding around a model with 7 ~ 1 in the 
redshift range z = [2 — 4], while the original analysis was 
based on simulations with 7 ~ 1.6. Furthermore, the flux 
power is computed using a Taylor expansion to second in- 
stead of first order. The parameters of the f iducial cosmolog- 
ical simulation are those of the B2 model in IViel fc Haehneltl 
(2006). As before the flux statistics have been corrected for 
box-size and resolution effects. We compute the derivatives 
required for the Taylor expansion by performing between 
four and six hydrodynamical simulations for every cosmo- 



logical and astrophysical parameter considered. In addition, 
we now allow for the effect of the reionization redshift, z Ie , 
and introduce this as an extra parameter. We interpolate be- 
tween two very different reionization histories with z IC ~ 15 
and 2 rG ~ 7. We also introduce two extra parameters de- 
scribing the redshift evolution of the thermal state of the 
IGM, the power-law index of the T and 7 relations at 2 > 3 
(a redshift range which is not probed by the PDF). 

The results for the power spectrum only analysis are 
summarised in the first two columns of Table [2] for a low 
redshift only sample, 2 = [2.2 — 3.6], and the full SDSS 
data set, 2 = [2.2 — 4.2]. We decided to perform a sepa- 
rate analysis whi c h omi ts the highest redshift bins following 
IViel fc Haehneltl (|2006h . who obtained a somewhat poorer 
fit for the high redshift (2 > 3.6) PS estimates. Despite the 
fact that the flux statistics were calculated by expanding 
around a reference model with very different thermal his- 
tory, in both instances the analysis still gives constraints on 
the cosmological pa rameters that are i n agr eement with the 
previous analysis of I Viel fc Haehneltl (|2006). This is rather 
reassuring. However, there are some aspects of the results 
that need scrutiny. First, we note that for the flux power 
only the temperature at mean density, To, is significantly 
higher than that preferred by the PDF (and higher than ex- 
pected for the photoionized IGM). Secondly, the value of as 
is now somewhat on the lower end of values allowed by the 
previous analysis of Lyman-a d ata and thus in better agree- 
ment with the CMB data (e.g. iKomatsu et all [2008). This 
is due to the degeneracy between as and 7 discussed in B08; 
allowing for 7 < 1 means that the flux power can now be 
fitted by a slightly lower as (but note that other parameters 
also have a a significant influence on the inferred as, most 
notably the mean flux level). 

Overall the results from the new flux PS analysis are 
consistent with those inferred from the PDF alone except 
for the value of To at 2 = 3, which is several a above that 
inferred from the best fit to the flux PDF. In the last col- 
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umn of Table [2] we show the constraints for a joint analysis 
of flux PDF and PS. Somewhat surprisingly the joint analy- 
sis prefers a larger value of as = 0.9 ±0.02 with rather small 
errors. We explicitly checked that this large value of as is 
related to the rather different To values that the PDF and 
PS favour. If we artificially remove the constraint of the tem- 
perature being simultaneously consistent with the somewhat 
discrepant temperatures favoured by the flux PDF and PS, 
the joint analysis gives as = 0.86 ± 0.03 with an improve- 
ment of A\ 2 = 12. A not yet accounted for systematic error 
in the measurement of the flux PDF and/or PS appears to 
be a possible explanation for this discrepancy. Alternatively 
the inconsistencies may suggest that a power-law T — p re- 
lation is a poor approximation to the thermal state of the 
IGM and a wider range of physically motivated relations 
should be considered in future simulations. 



5 DISCUSSION 

We have presented cosmological and astrophysical con- 
straints derived from the K07 Lyman-a flux PDF measured 
from a set of 18 high-resolution QSO spectra whose statisti- 
cal and systematic errors have been carefully estimated. The 
Lyman-a flux PDF on its own provides tight constraints on 
the thermal state of the IGM and on cosmological parame- 
ters describing the linear dark matter PS. The results have 
been obtained by fitting the flux PDF at three different red- 
shift bins in the range 2 < z < 3 and for three different 
flux ranges F = [0.1 — 0.8], F = [0 - 0.9], and F = [0 - 1]. 
There is good agreement between the analyses for the full 
flux range and the two restricted flux ranges and the re- 
sults are consistent with those derived in the simpler analy- 
sis made by B08. An inverted temperature-density relation 
is favoured at the ~ 3a level (at z ~ 3) if we consider the 
PDF for the full flux range, but the significance is reduced 
to 1 — 1.5a if we restrict the analysis to F = [0.1 — 0.8]. 
The constraints for other parameters are in agreement with 
th ose presented in the li terature (e.g. the SDSS flux power 
in iMcDpnald et all 2005 ) . We have also refined the method 
used by Viel fc Haeh nelt ( 2006J) and updated the constraints 
from the SDSS flux power spectrum. The constraints from 
the flux power spectrum are consistent with those from the 
flux PDF, with the exception that the flux power spectrum 
prefers a significantly larger temperature at mean density, 
To. A joint PS-PDF analysis gives a reasonable fit to the 
data but results in a larger as than an analysis of CMB 
data alone. This discrepancy appears to be related to the 
higher To that the flux power spectrum prefers. 

Recent simulations of photo-heating during 
He ii reionization indicate that an inverted T — p re- 
lation is very difficult to achieve within the standard 
model describing the thermal state of the IGM even 
if radiative transfer effects are taken into account, at 
least if Hen re i onizat i on is driven prima rily by QSOs 
l|McQuinn et all 120091 ; iBolton et all |2009a| ). It therefore 
appears difficult to reproduce the observed flux PDF with- 
out invoking a not yet identified source of IGM heating, 
or additional systematic errors which impact on the flux 
PDF. Similarly, the high To values preferred by the PS are 
very difficult to reconcile with constraints fro m the widths 
of thermally broadened absorption lines (e.g. ISchave et al.l 



(|200Ch ; [Ricotti et al.l (^POOD ) and our understanding of the 
heating of the photoionized IGM. 

As the results from the flux PDF and the CMB data 
agree very well, the rather high values of To preferred by 
the flux power spectrum suggests perhaps instead that we 
have identified a not yet accounted for systematic error in 
the SDSS flux power spectrum data. Independent analy- 
sis based on line statistics, on new data sets at medium 
and high-resolution and further progress in incorporating 
Hen reionization models into high resolution hydrodynam- 
ical simulations will hopefully allow us to further improve 
our understanding of the systematic uncertainties of Lyman- 
a forest data and resolve these small but statistically signif- 
icant inconsistencies. 
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